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1.  INTRODUCTION 


The  term  "inertial  navigation"  is  commonly  used  in  the  literature  to  describe  the 
process  of  obtaining  velocity  and  position  information  from  measurements  made  with 
inertial  instruments.  To  accomplish  this,  an  inertial  navigation  system  must  contain 
four  basic  conceptual  elements:  a  vector  accelerometer,  an  attitude  reference,  a 
computer,  and  a  clock  (Russell  1962).  Russell  describes  each  of  these  elements  as 
follows  (pp.  17-18).  A  vector  accelerometer  is  a  device  which  measures  the 
nongravitational  acceleration  of  the  platform  to  which  it  is  afixed.  Oftentimes  this 
vector  is  defined  by  three  orthogonal  components,  each  of  which  is  measured  by  a 
single-degree-of-freedom  (linear)  accelerometer.  The  attitude  reference  gives  the 
orientation  of  the  accelerometers  in  the  reference  frame  in  which  navigation  takes 
place.  Gyroscopes  have  most  often  been  employed  to  this  end.  The  computer  is  the 
conceptual  element  which  solves  the  acceleration  equation  by  calculating  gravitational 
acceleration,  by  performing  integrations,  and  by  making  coordinate  transformations  as 
required.  The  clock  is  needed  to  predict  the  gravitational  field  and  the  location  and 
orientation  of  moving  external  reference  frames  (e.g.,  the  earth  beneath  an  aircraft). 

It  has  been  shown  analytically  that  it  is  possible  to  obtain  the  information 
necessary  for  determining  both  linear  accelerations  and  angular  motions  using  only 
measurements  from  linear  accelerometers  (Schuler,  Grammatikos,  and  Fegley  1967). 
Efforts  are  underway  to  infer  the  gravity  vector  from  accelerometer  measurements. 
The  goal  is  to  develop  signal  processing  software  that  can  isolate  the  gravity-induced 
component  of  projectile  pitch  that  is  present  in  intra-atmospheric  trajectories. 
Because  high-g  miniature  accelerometers  are  common  but  gyros  are  not,  a  navigation 
system  which  uses  only  linear  accelerometers  for  inertial  measurements  would  be  an 
attractive  candidate  for  inclusion  in  artillery  projectiles.  Such  a  system  could  perhaps 
enable  such  things  as  autonomous  course  corrections  and/  or,  in  conjunction  with  a 
telemeter,  projectile  registration. 

As  a  start  in  investigating  such  possibilities,  equations  for  the  inertial  acceleration 
of  an  arbitrary  point  on  a  flight  body  have  been  derived  in  forms  appropriate  for 
inclusion  in  a  computerized  six-degree-of-freedom  (6dof)  trajectory  model  used  at  the 
U.  S.  Army  Research  Laboratory  (ARL).  This  trajectory  model,  called  CONTRAJ 
(Hathaway  and  Whyte  1991),  was  then  modified  to  allow  up  to  12  arbitrarily  located 
3-axis  accelerometers  on  the  projectile.  Using  appropriate  configurations  of 
accelerometer  locations  and  orientations  and  linear  combinations  of  accelerometer 
outputs,  the  components  of  projectile  linear  and  angular  motion  can  be  estimated. 
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Importantly,  it  was  found  that  the  configurations,  orientations,  and  combinations 
that  successfully  estimate  the  angular  motion  components  for  one  projectile  will 
sometimes  not  work  for  another  projectile  with  different  motion  characteristics. 
Specifically,  a  methodology  employed  by  Schuler,  Grammatikos,  and  Fegley  for 
estimating  the  components  of  angular  velocity  failed  due  to  numerical  problems  when 
implemented  on  the  M483A1  artillery  projectile.  An  alternative  approach  yielded 
accurate  estimates  demonstrating  that  accelerometer-only  inertial  measurement  units 
need  to  be  tailored  to  the  anticipated  flight  characteristics  of  the  body  on  v^ich  they 
will  be  installed.  Appropriate  accelerometer  configurations  and  computational 
algorithms  for  use  in  many  Army  projectiles  can  be  identified  using  the  computer 
model  developed  for  this  study. 


2.  LINEAR  ACCELERATION  OF  A  POINT 

Define  a  set  of  inertially  fixed  Cartesian  axes  with  origin  Oq  and  unit 

vectors  ^  q  o  ^  o  ^  second  set  of  Cartesian  axes  X  YZ  with  origin  O  and  unit 
vectors!^  7^^.  Not  only  may  O  move  relative  to  Oq  but  the  X  YZ  axes  may  rotate  as 
well.  Consider  a  point  H  moving  relative  to  both  sets  of  axes.  If  we  denote  H’s 
position  relative  to  in  the  X^Y^Z^  system  by  ff,  O ’s  position  relative  to  in  the 
XYZ  system  by  ^ ,  and  H ’s  position  relative  to  O  in  the  Z  FZ  system  by  It ,  then 

(1) 

and  it  can  be  shown  (e.g.,  Page  1952)  that  inertial  velocity  and  acceleration  oiH  are 
given  by: 

^  _  V  -  *  V 

—  (2) 

^+wxj^-l-wx(i;5x^) -f-2(u;  x^) (3) 

where  It  is  the  angular  velocity  of  the  XYZ  system  with  respect  to  the  XqYqZq 
system. 

The  velocity  is  the  sum  of  three  terms.  The  first  gives  the  velocity  of  the  origin  of 
the  moving  axes.  The  second  is  the  velocity  of  H  due  to  the  angular  velocity  w  of  the 
moving  axes.  The  third  ^  is  the  velocity  of  H  relative  to  O  as  measured  in  the 
moving  system. 
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The  acceleration  is  the  sum  of  five  terms.  The  first  is  the  the  acceleration  of  the 
origin  of  the  moving  axes.  The  second  is  the  linear  acceleration  due  to  the  angular 
acceleration  a;  of  the  moving  axes.  The  third  is  the  centripetal  acceleration  due  to  the 
angular  velocity  uJ  of  the  moving  axes.  The  fourth  is  the  Coriolis  acceleration  where 
^  is  the  velocity  of  H  relative  to  O  as  measured  in  the  moving  system.  The  fifth 
term  is  the  acceleration  of  H  relative  to  O  as  measured  in  the  moving  system. 

The  equations  of  motion  used  to  compute  projectile  and  rocket  trajectories  are 
usually  evaluated  in  either  of  two  right-handed  Cartesian  coordinate  systems.  The 
first,  called  the  body-fixed  system,  has  its  origin  at  the  center  of  mass  of  the  flight  body 
and  its  Z-axis  parallel  to  the  downrange  axis  of  symmetry.  The  V  and  Z  axes  are  then 
oriented  so  that  the  products  of  inertia  vanish.  The  three  axes  of  this  system  are  called 
the  principal  axes  of  inertia  of  the  body.  As  the  name  implies,  this  system  is  tied  to  the 
flight  body  and  translates,  orients,  and  rotates  with  that  body.  The  second  system, 
called  the  plane-fixed  system,  also  has  its  origin  at  the  center  of  mass  and  its  Z-axis 
along  the  axis  of  symmetry.  However,  the  Y  and  Z  axes  are  not  fixed  to  the  projectile. 
Rather,  the  T-axis  is  constrained  so  that  it  lies  in  the  horizontal  plane.  A  mutually 
orthogonal  Z-axis  then  completes  this  system.  Body-fixed  coordinates  are  usually 
employed  for  simulating  guided  flight  bodies  because  the  necessary  seekers,  sensors, 
and  control  mechanisms  are  fixed  to  the  bodies  and  most  easily  described  in  such  a 
system.  The  plane-fixed  system  eliminates  sensitivity  of  the  integration  to  projectile 
roll  rate  by  eliminating  the  affected  component  of  gravity.  For  spin-stabilized 
projectiles,  this  choice  of  coordinate  systems  can  result  in  significant  computer  savings 
by  not  requiring  the  extremely  small  integration  time  steps  necessary  to  avoid 
smearing  of  the  gravity  effect  over  the  roll  angle.  Though  it  is  possible  for  a 
sufficiently  large  flight  body  to  have  a  gravitational  gradient  within  the  body, 
CONTRAJ  assumes  a  constant  field  throughout  the  body.  In  CONTRAJ,  therefore, 
the  equation  of  linear  motion  is  d^=A  +~g  where  is  the  nongravitational 
acceleration  of  the  center  of  mass  (c.m.)  and^  is  the  gravitational  acceleration  of  the 
c.m.  A  is  also  called  by  some  authors  the  thrust  acceleration,  sensed  acceleration,  or 
specific  force. 

In  the  next  two  subsections,  the  derivations  of  the  body-fixed  and  plane-fixed 
system  versions  of  equations  (l)-(3)  are  summarized.  In  subsection  2.3,  the 
acceleration  components  as  computed  in  the  plane-fixed  system  will  be  rotated  into 
the  body-fixed  system  to  allow  for  comparison  with  the  analogous  estimates  made  for 
a  trajectory  computed  in  the  body-fixed  system.  Uninterested  readers  may  skip  these 
and  continue  with  section  3  without  loss  of  continuity. 
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2.1.  Accelerations  Computed  in  a  Body-Fixed  System. 


When  computing  a  trajectory  emplojing  the  body*fixed  coordinate  sy^em,  the 
vectors  in  equations  (l)-(3)  can  be  e^qjressed  as: 

=  AXo^’ 0 +Ay{)/ 0  0 

^  =axT  4-Ayy  +azI^ 
lit  =  ajcT  +  AyX  Az^ 

0^  =  pT  +  qX  +  rX  (4) 

oj  =  pT  +  qX  +  rX 
Xn  +Ay7^  +AZ? 

^  —  AX~f  +  AyX  +AzX 

where  p  and  p  are  the  angular  velocity  and  acceleration  of  the  j-k  axes  about  the  i 
axis,  q  and  q  are  the  velocity  and  acceleration  of  the  i-k  axes  about  the  j  axis^  and  r 
and  r  are  the  velocity  and  acceleration  of  the  i-j  axes  about  the  ^  axis. 

Making  these  substitutions  in  equation  (3),  forming  the  cross  products,  and 
collecting  terms  yields: 

ft  =  ^ -fT\{-q^-r')Ax  -\'{pq-r)Ay  -\'{pr  +  q)Az  +2{qAz-rAy)  +Ax] 

+r)Ax  +(-p^-r^Ay  -\-{qr~p)Az  +2(rAx-pAz)  +Ay]  (5) 

-\-X\(pr-q)Ax  +(qr  +p)Ay  +(-p^-q^)Az  +2(pAy  ~qAx)  -V-Az] 

The  components  of  velocity  of  the  c.m.  are  typically  represented  in  the  moving 
( X  y  Z  )  system  as  m  ,  v ,  and  w.  If  we  assume  that  H  is  fixed  on  the  body,  Ax ,  Ay ,  and 
Az  are  nonconstant  only  if  the  c.m.  is  not  fixed,  as  it  generally  would  not  be  if  fuel  or 
supplies  are  consumed  during  the  flight.  In  such  a  case,  it  is  useful  to  make  the 
following  substitutions: 


U  —  WQ+/Li 

v=Vq+t}  (6) 

w  =  Wq 
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where  Mq,  Vq,  and  Wq  are  the  components  of  velocity  of  the  point  that  was  the  initial 
c.m.  and  /x.r/,  and  ^  are  the  components  of  velocity  of  the  c.m.  due  to  the  c.m.’s 
changing  location  within  the  projectile. 

If  we  express  the  velocity  of  O  in  body-fixed  coordinates  as: 

&=(uQ+fi)T  +(Vq+77)/  +.(wq+0^  (7) 

it  follows  from  the  definition  of  (equation  [4])  that 

&=  (Uq~ax)T  +  (vQ-Ay)T  +  (8) 

and  after  differentiation: 

^  =  (Mq-Ax)7  +  (vQ-Ay)7^  +  (Wq-A2)^  +  uJx^  (9) 

Forming  the  cross  product  and  collecting  terms,  we  get: 

^  =  T|^(«q-ax) +^(>VQ-Az)-r(V(j-Ay)j 

+T[(Vo-A>')  +  r(tto-AJc)-p(wQ-Az)]  (10) 

+  r[(wQ-Az)  +p(vQ-Ay)-q{uQ-Ax)^ 

Finally,  substituting  (10)  into  (5),  the  body-fixed  components  of  H  are  given  by: 

ft.  =  [«q-ax -H  ^(wQ-Az)-r(vQ-Ay)]  (11a) 

+  +(pq-r)Ay  +(pr  +q)Az  +2(^Az -rAy) -l-Ax] 

=  [(VQ-Ay) +r(wQ-Ax)-;?(wQ-Az)]  (11b) 

+  [ipq  +r)AJc  -l-(-/?^-r^)Ay  +{qr~p)Az  -^2{rAx-pAz)  +Ay] 

^  +p(vQ-Ay)-q(uQ-Ax)^  (11c) 

+  \^r-q)AX  -h(^r+/7)Ay  +{-p^-q^)AZ  +2{pAy-qAX)  -hAz] 
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The  body-fixed  components  of  the  difference  in  the  accefetat^s  of  any  two 
points  and  are  independent  of  c.m.  motion  and  given  by: 

+(p^-r)(Ayj-Ay2)  -t-(pr +g)(AZ|-AZ2)  (12a) 

+/')(aj:i-A;c2)  +{-p^~r'){Ay^-£^y^  +{qr~p){d,z^-M^  (12b) 

=  (pr-q)iAx^-AX2)  +(qr-^p)(Ay^~Ay2)  +(-p^’q^)(AZy-AZ2)  (12c) 

2.2.  Accelerations  Computed  in  a  Plane-Fixed  System. 

The  plane-fixed  coordinate  system  with  its  origin  at  the  c.iii.  differs  from  the 
body-fixed  coordinate  system  with  its  origin  at  the  c.m.  by  the  angle  between  the 
plane-fixed  F-axis  and  the  body-fixed  F-axis.  This  angle  will  be  designated  by  A<^,  and 
the  axes  of  the  two  systems  will  be  defined  so  that  they  are  coincident  when  A<^  =  0. 
Denoting  the  distance  from  the  Z-axis  to  H  by  R  and  the  angle  between  the  F-axis 
and  the  projection  of  0-H  onto  the  F-Z  plane  by  the  vectors  in  equations  (l)-(3) 
can  be  expressed  as; 

it  =  AXqI  q+AYqJ  q+AZq^  q 

O  —  aX  i  +AFy  +  aZ  A’ 

it  =  Ax~t  +  R  cos(j>ffX  +  ^  ^bi(j>fflt 

uj  =  qj^  +  rlt  (13) 

OJ  =  rlt 

7t^  =  axT  -f-  |^/?cos(^^ -/?sin(^^(^J  7^  +  ^/?sin<^^  -l-/?cos(^^^^j  It 

it^  =  ax7  +  j^-/?^sin(^^(^^  +cos(^jy(^^j  -2Rsm(j)^^fj  +/?cos<^j7^ 

+  T/^fcosc^^^^  -sin(^jy(^^l  ^2Rcos(j)jj^jj  +R  sin<^^l  It 
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Making  these  substitutions  in  (3),  forming  the  cross  products,  and  collecting  terms  yields 

ft  =  ^  +T[Ax(-^^-r^)  +Rcos<j)^{2q(j)^-r)  +Rsm.(j)^{2r(j)^  +q) 

+  Rzos^fj{-lr)  •\-Rs,m<j)^{lq)  +Aj: 

+7^[AA:r  +Ax(2r)  +/?cos0^(-r^-0^)  +/?sin0^(/^-0^)  (14) 

+  Rsm.<f>j^{-2(j)jj)  +/?cos^^j 

+lt^Ax(-q)  +Ax(-2q)  +/?cos<^^(<yr +0^)  +Rsm(j)f^(-q^-^^) 

+  Rcos<j>jj{2(j)^)  +/?sin<^^j 

The  most  convenient  way  to  reconcile  the  body-fixed  {BF)  and  plane-fixed  {FP) 
representations  of  the  location  of  the  point  H  is  to  synchronize  the  two  systems  initial 
conditions.  With  the  coincidence  of  the  two  systems  (A0  =  O),  the  respective 
representations  are  component -wise  equal.  With  AXgp=AXpp,  Aygp -Rcos<j)Q, 
AZgp  =Rsm(j>Q  initially  and  denoting  the  cumulative  rotation  of  the  flight  body  by  (j), 
the  following  relations  hold: 

^  ^0 
(j>^  =  <(>+(j>Q 

Rcos(j)Q  =  Ay  +Az^q  (15) 

Rsm<f)Q  -  Az  -Ay^Q 

••  ..  -  2  '  *  ■■ 

Rcos(f)Q  =  Ay  -  Ay(f)Q  +2az<^o 

Rsm(j>Q  =  AZ  -  az<I>q  -  2Ay<j>Q  -  Aycj)^ 

With  these  substitutions  into  (14)  and  sufficient  diligence,  the  components  of  ^are 


then: 


ft.  =  o'.  +Axi~q^-r^)  +AX 


(16a) 


+  Ay[cos<^(2^(^r)  -l-sin0(2r<^-l-^)]  -l-Ay(2^sin(^-2rcos<^) 
-l-Az[cos<^(2r(^-l-g)-sin(^(2^(^r)]  -l-Az(2<?cos<^  +2rsin<^) 
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ttj  +AJcr +AJc2r 


(m) 


-hAy\cos<l)(~r^~^)  +Ay(-2^ia^)  +Aycos^ 

+ Az\cos(j)(rg  -<l>)  - sin0(  -r^ +  Az{  *2^0$^)  +  AZ ( -sin<^) 

^  ^  +AJc(-^)  +ajc(-2^)  (i6c) 

+  Ay[cos0(r(9  +^)  +aj(2^^^)  -J-Ay  sin^ 

+Az[cos<^-<7^-(^^)-sm4^/^  +<^)]  +Az(-2^sjii^)  +AZ(cos^) 

Finally,  letting p  =<j>  and  p  =(j>,  the  plane-fked  components  of  the  acceief4tion  of  B 
in  the  most  widely  used  notation  are  given  by: 

^  ^  +A:r(-^^-r^)  +Ax  (17a) 


(17a) 


+Ay[cos^2^ -r)  +sin^2/p  -f^)]  4'Ay(2fsih^-2?^e0®^) 

4-Az[cos<^2/p  +^)-sin<^2^-r)]  +Az(2^cos^  +  2rsin^) 
itj  =  d^j  +A^r  +AJr2r  (17b) 

+Ay[cos<f)(-r^-p^)  +sin(^A^-p)]  +Ay(“2ps!n<^)  +Aycos0 

+Az[cosj^/9 -p)-sin(^-r^ -p^  4-Az('2pcos^)  +Az(*SHi^) 

^  ^  +Aa:(-^)  +a;c(-2^)  (17c) 

+  Ay[cos<^r<?  +p)  +sin<^>^^-p^)]  +Ay(2pcos^)  +Ay  sin^ 
+Az[cos(^-^^-p^)-sin<^7^  +p)]  +Az(-2psin^)  *f  Az(cos^) 

2.3.  Transformation  of  Acceleration  Components  From  Plane*roed 
to  Body-Fixed  System. 


Thus  far  the  location  of  H  has  been  reconciled  in  b^ 
and  expressions  have  been  derived  for  the  coinponeiits  of 


at  M  HI  both 


systems  (equations  [11]  and  [17]).  It  remains  to  rotate  the  components  into  a  common 
system.  Since  strapped-down  accelerometers’  measurements  will  be  at  fixed 
orientations  with  respect  to  the  flight  bodies,  the  plane-fixed  components  will  be 
rotated  into  the  body-fixed  system  using  the  following: 


(ft.  1 

AO  0  > 

(ff.  ) 

'SF 

tFP 

— 

0  cos<t>  sm<f> 

% 

(18) 

ff 

<.0  -sintj)  cos<j)^ 

ffkFp' 

Performing  the  matrix  multiplication  and  equating  components  gives: 


ff.  -  H.  cosS+Hj^  sm<j) 

JBF  JfP  ^FP 

ft.  = -ff;  sm(j)-\-ff.  cos<j> 

^BF  JFP  ^FP 


(19a) 

(19b) 

(19c) 


After  substituting  (17)  into  (19),  collecting  terms,  and  employing  several  trigonometric 
identities,  the  acceleration  components  in  the  body-fixed  system  in  terms  of  the 
angular  motion  components  in  plane-fixed  coordinates  are  given  by; 


ff.  -ff  +Ax(-^^-r^)  +AX 


(20a) 


JBF 


+  Ay[cos0(2qp -r)  +sin0(2;p  +q)]  +Ay(2qsin^-2rcos(^) 

+Az[cos0(2;p  -l-q)-sin<^2qp-r)]  +Az(2qcos<^  +  2rsin(;i) 
^.cos<^  +  ^sin<^ 

+  Ajc(rcos0-^sin^)  +Ax(2rcos(^-2^sin(^) 


(20b) 


+Ay 


-P 


r  2  .  2> 

r  -\-q 

r  2  2'k 
r  -q 

.  2  . 

.  2  > 

cos2<^+^rsin2<^ 


+Ay 


+  AZ 


-P  + 


r  2  2^ 

r  -q  ' 


sin2(^  +  ^rcos2(^ 


+Az(-2/7) 
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(20c) 


it,.  = 


'•BF 


d'j  ( sincj))  +  icos(f)) 

+ AX  ( -r  sin  <j)-qcos(l>)  +Axi  -2rsm  (j>-2qc<^S(l)) 


H-Ay 


+  AZ 


P  + 


2 

-P  - 


2  2>, 

r  -q 


sin2<^  +  ^rcos2^i 


+Ay(2/7) 


2  ,  2\ 
r  +q 


t  2 
r  -q 


cos2(f)~qrsm2m  -hAZ 


The  body-fixed  components  of  the  difference  in  the  accelerattons  of  any  two 
points  and  in  terms  of  the  angular  motion  components  in  plane-fixed 
coordinates  are  given  by: 


‘bf  ^bf 


ft.  -ft.  =  (-^^-r^](AXj-AX2) 

[cos<^2^-r)  -hsin^lip  +^)]  (Ay  Ay 2) 

\cos(j){2rp  +q)-sm(f>i2qp -f^  {AZi~AZ2) 

ft.  -ft  =  [  rcos^  -  ^  suk;^]  (  AX  j  -  AX2) 

j  1  j  ^ 


(21a) 


-h 


-I- 


(21b) 


^BF  -^BF 


+ 


-1- 


2 

-p  - 


-p  4- 


r  2  .  ^2>| 
r  +q 

r  2  2\ 
r  -q 

.  2  J 

.  2  . 

2  2^ 

r  -q 


cos2<^ +^rsm2^l 


(Ay  Ay  2) 


I  2 


sin2<^  +  qrcosTxm  (Az  1  -  AZ2) 


^  =  [-rsin(^-^cos<^](AXj-AX2) 


(21c) 


^BF  ■‘'BF 


-h 


-h 


r  / 


V  + 


2  2\ 
r  -q 


I  2  ) 


sin2<^  +  qrzos  2(m  (Ay  ^  -  Ay2) 


2 

-p  - 


(  2 


f  2  2l 

r  -g 


cos2<^-^/'sin2^i 


(AZj-AZ2) 
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3.  SIMULATION  RESULTS 


Though  strapped-down  linear  accelerometers  could  be  installed  anywhere  on  a 
projectile  and  oriented  arbitrarily,  the  algebra  required  to  derive  rotational  motion 
components  from  the  accelerometer  measurements  is  greatly  simplified  if  the 
measurement  axes  of  the  accelerometers  are  oriented  parallel  to  the  principal  axes  of 
inertia  of  the  projectile  (Figure  1).  An  accelerometer  so  positioned  would  be  directly 
measuring  one  of  the  body-fixed  components  of  nongravitational  inertial  acceleration 
at  that  location.  Thus,  an  accelerometer  at  ax,  Ay,  Az  from  the  projectile  c.m. 
oriented  in  the  47^  direction  would  be  measuring  ^  -  ^  and  similarly  for  +7^  and  -\-lt 
orientations.  When  these  components  are  expressed  in  terms  of  the  linear  and 
angular  motion  components  in  body-fixed  coordinates  (equation  [11])  and  the 
difference  between  the  accelerations  at  two  different  locations  is  similarly  expressed 
(equation  [12]),  methods  for  deriving  the  components  of  angular  velocity  and 
acceleration  can  be  readily  demonstrated. 

tt.  =  [^M(j-Ax +^(W(j-Az)-r(vQ-Ay)]  (11a) 

+  [(-^^-r^)Ax  +{pq-r)Ay  +(pr-l-<7)AZ  +2(^Az-rAy)  +Axj 
ttj  =  [^(Vg-Ay) +/-(Mo-ax)-;7(Wq-az)]  (11b) 

+  [(P^  +r)Ax  -l-(-p^-r^)Ay  +(^r-/7)AZ  +2(rAx-/7Az)  +Ay] 

^  [(Wq-Az)  +p(VQ-Ay)-g(MQ-AX)j  (11c) 

+  [(pr-^)Ax  +  {qr  +/7)Ay  -h(-p^-<?^)Az  +2(/?Ay-^Ax)  +Az] 

=  (-^^-r^)(AXj-AX2)+(p^-r)(Ayi-Ay2)+(p^+^)(AZi-AZ2)  (12a) 

=  (p^ +r)(AXj-AX2)+(-p^-r^)(Ayi-Ay2)+(^r-/7)(AZi-AZ2)  (12b) 

=  (pr-^)(AXj-AX2)+(^r-hp)(Ayi-Ay2)+(-p^-^^)(AZi-AZ2)  (12c) 

where  p  and  p  are  the  angular  velocity  and  acceleration  of  the  j~k  axes  about  the  i 
axis,  q  and  q  are  the  velocity  and  acceleration  of  the  i~k  axes  about  the  j  axis,  and  r 
and  r  are  the  velocity  and  acceleration  of  the  i-j  axes  about  the  k  axis. 
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Schuler,  Grammatikos,  and  Fegley  postulated  a  configuration  of  six 
accelerometers  positioned  as  follows: 


=AX+7^,A>',A^ 
=  AX  +7^,Ay,Az 
=  Ax,Ay  +7^,Az 
=  Ax,A>'  +7^»Az 
=  Ax,Ay,Az  +7^ 
//y  =  AX,  Ay,  Az +7y. 


aligned  with  the  i  axis 
aligned  with  theT  axis 
aligned  with  they^  axis 
aligned  with  the 7^  axis 
aligned  with  the  k  axis 
aligned  with  the  it  axis 


Denoting  the  corresponding  sensed  accelerations  by  and  A ,  they 

obtained: 


or 


similarly 


and  ultimately 


2  .  2. 


{q  +r)  = 


D  a 


x-% 


s  C 


A.  ~A^ 

a  c 


7.-7/ 


=  c. 


p  = 


q  = 


-Cj  +C2  +  C3 


Ci“C2  +  C3 


r  = 


(22a) 


(22b) 


(22c) 


(23a) 

(23b) 

(23c) 
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They  point  out  the  sign  difficulty  in  obtaining  p ,  q,  and  r  and  suggest  that  this 
"difficulty  can  be  resolved  through  the  use  of  auxiliary  devices  that  may  be  less 
accurate  and  less  costly  than  accelerometers,"  Implementation  of  this  configuration 
and  computational  scheme  in  a  simulated  flight  of  an  M483A1  artillery  projectile 
revealed  an  additional  difficulty  with  this  methodology  for  such  projectiles. 

3.1.  Case  1:  M483A1  Artillery  Projectile. 

The  M483A1  is  a  155-mm-diameter,  spin-stabilized  projectile  that  is  launched 
from  a  rifled  tube  with  a  1/20  twist.  When  fired  from  an  M109  howitzer  with  a  7w 
propelling  charge,  the  M483A1  has  a  launch  velocity  of  539  m/ s  and  an  initial  spin  rate 
of  ==174  rps.  At  45^  launcher  elevation,  under  nominal  conditions,  this  projectile  has  a 
range  of  =14  km  and  a  flight  time  of  =60  s.  At  impact,  the  spin  rate  is  =125  rps. 
Thus  p,  whose  units  are  radians/ s,  ranges  from  =1100  rad/ s  to  =785  rad/ s 
(Figure  2a).  At  the  same  time,  q  and  r  range  in  magnitude  from  near  zero  at  launch 
to  =0.6  rad/  s  at  impact  (Figures  2b  and  2c). 

Though  the  Schuler,  Grammatikos,  and  Fegley  method  for  finding p,q,  and  r  is 
analytically  correct,  it  sometimes  failed  to  produce  accurate  values  of  q  and  r  for 
simulated  flights  of  this  projectile  due  to  numerical  difficulties.  Whereas  p  is  on  the 
order  of  magnitude  of  10  ,  ^  and  r  are  on  the  order  of  magnitude  of  10  or  less.  Thus 
there  is  an  at  least  four  orders  of  magnitude  difference  between  p  and  q  and  p  and  r 
and  an  eight  orders  of  magnitude  difference  between  their  squares.  Significant  figure 

limitations  then  result  in  estimation  errors  for  the  values  of  q  and  r  due  to  the  Schuler, 

2  2 

Grammatikos,  and  Fegley  methodology’s  reliance  on  accurately  determining  (p  -\-q  ) 
and  (p^  +r^. 

An  M483A1  trajectory  with  launch  conditions  as  described  above  was  computed 
with  the  known  values  of  p,  q,  and  r  and  their  Schuler,  Grammatikos,  and  Fegley 
estimates  output  to  a  file  at  0.015-s  intervals  throughout  the  flight.  Additionally, 
estimates  of  q  and  r  generated  using  a  different  combination  of  accelerometer 
positions  and  a  different  computational  algorithm  were  output.  Within  the 
CONTRAJ  model  as  implemented  on  our  branch  computer,  64  bits  are  used  to 
represent  reals,  52  bits  of  which  are  used  for  the  mantissa,  11  for  the  exponent,  and  1 
for  the  sign.  This  gives  15  significant  figures  to  the  known  values  ofp,  and  r  and 
the  accelerometer  outputs. 

Figure  3  gives  the  percentage  of  the  samples  for  which  the  errors  in  the  estimates 
of  q  and  r  exceeded  each  of  four  different  levels  as  a  function  of  the  number  of 
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significant  figures  carried  in  the  computations  for  both  methodok)gies.  Figure  3a  gives 
the  percentages  of  time  that  the  magnitudes  of  the  error  hi  the  estanate  of  q  and  r 
exceeded  5%  of  the  known  magnitudes  of  q  and  r.  Figure  3b  gives  the  percentages 
for  which  the  error  magnitudes  exceeded  25% ,  3c '  50% ,  and  3d  >  100% .  The  upper 
set  of  curves  in  each  plot  is  for  the  errors  in  the  Schuler,  Grammatikos,  and  Fegley 
estimates  of  q  and  r.  Only  at  14  ai^  15  significant  figures  are  these  estimates 
approximately  errorless.  The  curves  for  the  alternative  estimates  of  q  and  r  are 
effectively  colinear  and  are  seen  as  the  lower  curve  in  each  plot.  These  estimates  are 
nonzero  only  for  two  significant  figure  computations.  In  this  case,  the  estimation  error 
is  primarily  due  to  round-off  as  the  estimates  are  always  being  compared  to  the  15 
significant  figure  known  values  of  q  and  r. 

Though  this  alternative  method  of  estimating  angular  velocky  components  has 
the  advantage  of  requiring  less  significant  figures  in  its  computations  (and 
accelerometer  data),  it  has  the  disadvantage  of  requiring  the  output  of  10 
accelerometers  rather  than  6.  Their  configuration  is: 


=Ax,Ay +7^,Az 
=A;r,Ay +7^,Az 
=  Ax,Ay,Az  +7^ 
=  Ax,Ay,Az  +7^ 
=  Ax+7^,Ay,Az 
=  AX  -1-7^,  Ay,  Az 
=  AX  +7^,Ay,Az 

=  Ax+7^,Ay,Az 
H.  =  Ax,Ay  +7.,  AZ 
Hj  =  Ax,Ay  +7y,Az 

If  the  corresponding  sensed  accelerations 
and  i4y,  then: 

.  a  I 

pq-r  =  - 

la-lb 

.  c  a 

pr+q  =  - 

Ic-ld 

. 

pq+r  =  - 

le-lf 


aligned  with  the  i  axis 
aligned  with  theT  axis 
aligned  with  theT  axis 
aligned  with  theT  axis 
aligned  with  theT  ^xis 
aligned  with  the  j  axis 
aligned  with  the  T  axis 
aligned  with  the  k  axis 
aligned  with  theT  axis 
aligned  with  theT  axis 


=  Cj  (24a) 

=  C2  (24b) 

=  C,  (24c) 


are  denoted  A^,Ag,Af,Ag. 
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pr-q 


(24d) 


A, -A. 

p^+?  =  -  =  Cj  (24e) 

'ij-'ii 

2  2 
Recalling  that  p  is  eight  or  more  orders  of  magnitude  greater  than  r  ,  the 

approximation ^p^  +r^  is  made.  The  estimates  of  the  angular  velocity  components 

are  then: 


p  =  (25a) 

C1+C3 

q  =  -  (25b) 

C2+C4 

r  =  -  (25c) 

Also  readily  available  are  two  analytically  exact  measures  of  angular  acceleration 
components: 


r  = 


^3-^1 

2 

^2-^4 

2 


(26a) 

(26b) 


3.2.  Case  2:  2.75-Inch  MK66  Rocket. 


The  preceding  computations  were  also  made  for  a  representative  trajectory  of 
another  projectile  with  very  different  kinematics,  the  2.75-in  MK66  rocket.  For  this 
trajectory,  the  motor  was  modeled  as  burning  for  1,12  s,  at  the  end  of  which  the 
projectile  speed  was  =720 m/s.  When  launched  at  45^^  elevation  under  nominal 
conditions,  the  time  of  flight  was  =50.6  s  and  the  range  to  impact  was  =8050  m.  The 
known  values  of  the  angular  velocity  components  /?,  ^ ,  and  r  are  seen  in  Figures  4a, 
4b,  and  4c  respectively.  For  this  trajectory,  the  differences  between  the  order  of 
magnitude  of  p^  and  those  ofq^  and  are  much  less  than  those  seen  for  the  M483A1. 
Plotting  the  percentages  for  which  the  estimation  errors  exceeded  the  four  levels  as 
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before  (Figure  5),  the  curves  for  errors  in  q  and  r  are  essentially  colinear  for  both 
methodologies  for  this  projectile.  Also,  the  character  of  the  curves  is  seen  to  be  like 
those  for  the  M483A1  albeit  at  different  (lower)  values  for  significant  figures.  The 
Schuler,  Grammatikos,  and  Fegley  methodology  (upper  curves)  is  outperformed  by 
the  alternative  methodology  at  all  levels  for  computations  maintaMimg  tess  than  eight 
significant  figures.  The  lower  curves,  like  their  counterparts  for  tlje  M483AI,  are 
constant  valued  except  for  the  single  instance  of  two  significant  figtares  and  estimation 
errors  exceeding  5% .  It  is  not  apparent  because  of  the  scale,  but  the  oonitaut  value  is 
not  -zero  as  it  was  for  the  M483A1.  It  is  0.653  % .  This  k  due  to  the  Renter  error  in 
the  approximation  of  p  +r  .  This  approximation  can  be  ehminated  if  12 
accelerometers  are  employed  rather  than  10.  With  through  through 

and  Cj  through  as  before,  and 


H.  =  Ax,Ay,Az  +7- 


Hj  =  AX,  Ay,  Az  +7y 
=  Ax,Ay  +7;t,Az 
=  AX,  Ay  +7p  Az 

. 

qr-p  =  - 


7,-7; 


aligned  with  the  j  axis 

aligned  with  the  X  axis 

aligned  with  the  Ic  axis 
aligned  with  the  it  axis 


qr  +p 


— -  =  c, 

7* -7, 


(27a) 


(27b) 


The  angular  velocity  components  are  then  given  by; 


P  = 


^  = 


r  - 


(Cj  +C^)(C2  +  C^)^ 


1/2 


2{C^  +  C^) 


Cj+C3 


2p 

^2  ^4 


2p 


and  the  angular  acceleration  components  by: 


^6-^5 


P  = 


(28a) 


(28b) 

(28c) 


(29a) 
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Q 


(29b) 


C  -C 
2 

C  -C 

r  =  -  (29c) 

2 

In  this  methodology,  the  sign  difficulty  analytically  remains  only  for  p .  In  practice,  this 
difficulty  appears  only  for  projectiles  whose  spin  direction  is  either  unknown  or  subject 
to  reversal  during  flight.  In  these  instances,  an  auxiliary  sensor  could  be  used  to 
remove  the  ambiguity  as  suggested  by  Schuler,  Grammatikos,  and  Fegley. 

4.  SUMMARY 

Equations  for  the  inertial  acceleration  of  an  arbitrary  point  on  a  flight  body  have 
been  derived  in  the  two  coordinate  systems  most  commonly  employed  in 
computerized  trajectory  models,  the  body-fixed  and  plane-fixed  systems.  The 
equations  were  constructed  in  such  a  form  that  the  location  of  the  point  was 
reconciled  in  both  systems.  The  acceleration  components  in  the  plane-fixed  system 
were  projected  onto  the  principal  inertial  axes  of  the  flight  body  for  comparison  with 
the  components  as  computed  in  the  body-fixed  system.  Each  of  the  acceleration 
components  given  by  these  equations  represent  the  sum  of  the  acceleration  that  will 
be  measured  by  a  perfect  accelerometer  located  at  that  point  on  the  flight  body  whose 
measurement  axis  is  oriented  parallel  to  the  respective  principal  axis  of  inertia  and  the 
projection  of  the  gravity  vector  onto  that  axis.  The  equations  for  these  components 
([11]  and  [17])  further  show  that  these  accelerations  are  given  by  various  combinations 
of  the  translational  and  angular  velocities  and  accelerations  that  constitute  the  body’s 
kinematics.  These  equations  were  included  in  the  CONTRAJ  6dof  trajectory  model. 
By  appropriately  locating  and  orienting  a  collection  of  accelerometers  and  forming 
linear  combinations  of  their  outputs,  solutions  for  the  velocity  and  acceleration 
components  can  be  derived. 

Implementation  of  a  configuration  of  accelerometers  and  a  computational 
algorithm  published  in  a  paper  by  Schuler,  Grammatikos,  and  Fegley  revealed  a 
further  consideration  in  the  design  of  accelerometer-only  inertial  measurement 
devices.  Though  the  Schuler,  Grammatikos,  and  Fegley  method  for  finding  angular 
velocity  components  is  analytically  correct,  it  sometimes  failed  to  produce  accurate 
estimates  of  these  velocities  in  a  simulated  flight  of  an  M483A1  projectile  due  to 
numerical  difficulties.  The  successful  estimation  of  these  velocities  employing  an 
alternative  approach  demonstrated  both  the  necessity  of  application  driven  design  of 
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such  systems  and  the  utility  of  the  accelerometer  augmented  CONTRAJ  model  in 
such  efforts. 

Realistic  models  of  accelerometers  need  to  be  added  to  CONTRAf  as  a  next 
step  in  investigating  the  possibility  of  practical  accelerometer-only  inertial  systems. 
How  bias,  noise,  and  other  corruptions  of  the  ‘perfect’  acceleratwn  measurements 
assumed  in  this  study  impact  such  systems  will  be  the  subject  of  ftiture  efforts. 
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ACCELEROMETERS 


Figure  1 .  Accelerometer  locations  and  orientations.. 


Percentage  Percentage 


c.  Estimate  Error  Exceeds  50% 


d.  Estimate  Error  Exceeds  100% 


Figure  3.  Angular  motion  estimation  errors  -  M483A1  projectile. 


rad./ sec  q  —  rad/sec  p  —  rad/sec 


Time  (sec) 


b.  Angular  Velocity  (q)  of  i-k  Axes  About  j  Axis 


Time  (sec) 

c.  Angular  Velocity  (r)  of  i-j  Axes  About  k  Axis 

Figure  4.  Bociy-fixed  angular  velocity  components  -  2.75-in  MKS6  rocket. 
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percentage  Percentage 


a.  Estimate  Error  Exceeds  5% 


b.  Estimate  Error  Exceeds  25% 


c.  Estimate  Error  Exceeds  50%  d.  Estimate  Error  Exceeds  100% 

Figure  5.  Angular  motion  estimation  errors  -  2.75-in  MK66  rocket. 
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